Hole-doping effects on a frustrated spin ladder 
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Hole-doping effects are investigated on the t-J ladder model with the linked-tetrahedra structure. 
We discuss how a metal-insulator transition occurs upon hole doping with particular emphasis on 
the effects of geometrical frustration. By computing the electron density and the spin correlation 
function by the density matrix renormalization group, we show that strong frustration triggers a 
first-order transition to a metallic phase, when holes are doped into the plaquette-singlet phase. 
By examining spin excitations in a metallic case in detail, we discuss whether the spin-gap phase 
persists upon hole doping according to the strength of frustration. It is further shown that the 
lowest excited state in a spin-gap metallic phase can be described in two independent quasiparticles. 



I. INTRODUCTION 



II. FRUSTRATED LADDER MODEL 
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Quantum spin systems with geometrical frustration 
have attracted considerable attention recently. There 
are a number of typical rnaterials, such as the spin- 
gap compQund SrCu2(B03)2a with the orthogonal-dimer 
structureJjIj'Ll the pyrochlore compounds with the lat- 
tice iDLjCQmerpaharing tetrahedral network of magnetic 
ions,ffl|ffle etc. 

Besides frustrated spin systems studied intensively 
so far, metallic systems with geometrical frustration 
have attracted much interest. For instance, it has 
been claimed that frustration due to a tetrahedral net- 
work of V ions may be important to understand heavy- 
fermion .bshavior of the metallic pyrochlore compound 
LiV204.E3 Also, i-another prototypical pyrochlore com- 
pound Y(Sc)Mn2ll^ shows spin-liquid behavior in a metal- 
lic phase. 

Stimulated by these experimental findings, we study 
here the hole-doping effects on fully frustrated spin sys- 
tems. To this end, we consider a simplified model in 
one dimension (ID), which still possesses the essence of 
strong frustration in the above-mentioned systems. To be 
precise, we employ the i-/ ladder model with the orthog- 
onal dimer structure, which also has linked-tetrahedra 
as a key structure. The advantage to employ the ID 
model is that this simplification allow us to accurately 
calculate physical quantities with use of such a reliable 
numerical method, as the density matrix renormalization 
group (DMRG),ll3 and thus precisely study the effects of 
frustration on hole-doped systems. 

This paper is organized as follows. In Sec. II, we in- 
troduce the frustrated t-J ladder model, and then in Sec. 
Ill study a metal-insulator transition upon hole doping. 
By computing spin correlation functions, we clarify the 
origin of the first-order phase transition found there. In 
Sec. IV, spin excitations in a metallic phase are dis- 
cussed in detail. It is clarified to what extent a spin-gap 
metallic phase persists around the insulating phase. In 
particular, we show that a spin triplet excitation in the 
spin-gap phase can be described in terms of two indepen- 
dent quasiparticles. We obtain the phase diagram in Sec. 
V, and summary and discussions are given in Sec. VI. 



The model we treat here is a frustrated electron lad- 
der which is an extension of the spin-1/2 Heisenberg 
model witli the orthogonal-dimer structure introduced by 
Gelfand.lla The model is illustrated schematically in Fig. 
|l|. It has two competing antiferromagnetic interactions 
Ji and J2, giving rise to strong geometrical frustration. 
In a hole-doped case, we use the t-J model where electron 
hoppingis assumed to have two distinct values, as shown 
in Fig. hi. The Hamiltonian is thus given by 



H 



E 



7n^nj ,(1) 



where the summation is taken over nearest neighbor 
bonds («,j) indicated by the thick and thin lines in Fig. 
|l|. Here, Cia is the annihilation operator of an electron 
with spin a at the i-th site and rii (Si) is the electron 
number (spin-1/2) operator. As mentioned above, hop- 
ping ti (exchange coupling J^) takes either ti or t2 (Ji 
or J2). Note that doubly occupied state at each site are 
implicitly forbidden in the t-J model. 
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FIG. 1: Frustrated ladder system, which is often referred to 
as a linked-tetrahedra chain. The thick and thin lines corre- 
spond to the parameters Ji (ti) and J2 (t2), respectively. 



We note here that the above modelll^MO'Ea has sev- 
eral remarkable properties in addition to those for or- 
dinary ID frustrated spin systems: (i) it is a ID 
analogue of the 2|D-| orthogonal-dimer model (Shastry- 
Sutherland modeffl) relevant for SrCu2 (603)2, for 
which the ground state is exactly given by a direct prod- 
uct of decoupled dimers; (ii) it is also regarded as a 
ID variant of the pyrochlore system (so-called linked- 
tetrahedra chain), which has frustrated tetrahedra as a 
key structure JI3 A remarkable point related to (ii) is that, 
for noninteracting electrons on this lattice, the model has 
a dispersionless fiat-band mode. 



e{k) = 



'tl — 4^2 cos(fc). 



(2) 



which reflects special geometry of the system, as is the 
case for the 3D pyrochlore lattice. Therefore, using this 
simplified model we can also explore the role of a flat- 
band mode, which is important to understand properties 
of the pyrochlore system near half filling. 



III. METAL-INSULATOR TRANSITION 

Let us start our discussions with the ground-state 
properties of the t-J model in a metallic phase. We par- 
ticularly focus on the nature of a metal-insulator tran- 
sition induced by hole doping. Here, we use a variant 
of the infinite DMRG method, so-callednproduct-wave- 
function renormalization group mcthodli^, to compute 
several quantities in a bulk limit. 

In the undoped spin model, the system is either in 
the exact dimer-singlet phase (J2/J1 0.71) or the 
plaquette-singlet phase {J2/J1 > 0.71).t3 We thus dis- 
cuss the hole-doping effects on these phases separately. 
In the following analyses, we set t2/ti = J2/J1 = 7, and 
mainly show the results of Ji/<i = 0.25 {ti = 1.0) for 
simplicity. 



A. dimer phase 

We begin with the dimer phase, which is stabilized at 
half filling for the ratio of the exchange couplings, 7 < 
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FIG. 2: (a) Electron density n and (b) short-range spin 
correlations Sa as a function of the chemical potential /i for 
7 = 0.3 at Ji/ti = 0.25. The system is in a dimer-singlet 
insulating phase at half filling n = 1. 



0.71. The calculated electron density n is shown in Fig. 
^(a) as a function of the chemical potential /j. 

It is seen that hole-doping smoothly drives the system 
to a metallic phase. At quarter filling {n = 1/2), there 
appears a plateau in n, which is due to formation of a 
CDW insulating state. It is instructive to compare this 
smooth transition with first-order transitions found in the 
magnetization curve at half filling in the dimer phase:E£l 
the magnetization exhibits jumps among the plateaus at 
zero, half, and full magnetization. The discontinuity in 
the magnetization reflect the fact that a triplet excited 
state in the dimer phase is completely localized due to the 
orthogonal-dimer structure. This is not the case for hole 
doping: doped holes are indeed heavy but still mobile, 
allowing a smooth transition to the metallic phase. 

We next compute spin correlation functions, which are 
defined in the rung (_L) and chain (||) directions as, 

S± — {SIS2), 



Shown in Fig. |2|(b) are the computed spin correlation 
functions. As the system approaches half filling from a 
metallic side, along the rung rapidly decreases down 
to —1/4 whereas along the chain becomes almost zero. 
Recall that the values of 5^ — —1/4 and S\\ — are ex- 
pected for a system with isolated dimers. Therefore, the 
above results of correlation functions imply that the sys- 
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first-order phase transition observed above. 



We note here that such a first-order metal-insulator 
transition only occurs when holes are doped in the pla- 
quette phase, but not in the dimer phase. It is also 
instructive to compare this with the results known for 
the magnetization curve at half filling, where first-order 
phase|-transitions occur only in the dimer phase (7 < 
0.71).t3 These apparently different aspects in phase tran- 
sitions come from the fact that hole doping induces a 
competition between the dimer- and plaquette-singlet 
states, whereas a magnetic field induces a competition 
between the dimer- (or plaquette-) singlet state and the 
triplet state. 
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FIG. 3: (a) Electron density n and (b) short-range spin corre- 
lations Sa as a function of the chemical potential /i for 7 — 0.9 
and Ji/ti = 0.25. This corresponds to hole doping in the 
plaquette phase. The dotted line in (a) indicates a jump as- 
sociated with a first-order transition. 



tem near half filling is regarded as a resonating metallic 
state of decoupled dimers. 



B. plaquette phase 

When holes are doped into the plaquette-singlet phase, 
the nature of a metal-insulator transition is completely 
changed. Shown in Fig. ||(a) is the electron density as a 
function of the chemical potential. It is remarkable that 
the system is driven to a metallic phase via a first- order 
phase transition, which is characterized by a jump in the 
electron density. 

We can see the origin of the first-order metal-insulator 
transition by checking short-range spin correlation func- 
tions in the rung as well as leg directions, shown in Fig. 
|(b). At half filling, 5^ > and < 0, being consis- 
tent with the plaquette-singlet ground state. However, 
once the system is driven to a metallic phase via the 
first-order transition, both of S± and become neg- 
ative. Furthermore, when the system approaches half- 
filling from a metallic side, S± decreases with negative 
values while S\\ gets very small. Namely, in a metallic 
phase close to half filling dimer correlations along the 
rung are largely enhanced. Hence, in the vicinity of half 
filling, this dimer-dominant state strongly competes with 
the plaquette-singlet state energetically, giving rise to the 




FIG. 4: Spin gap As and binding energy As of two holes 
{Nh = 2) for Ji/ti — 0.25. Inset shows a finite-size scaling 
analysis of 2 x L systems for 7 = 0.6. 
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''^ and spin density (S'f ) for two holes 
: 0.3 for (a) and (b) while 7 = 0.6 



FIG. 5: Hole density (n:. 
on the 2 X 64 system: 7 
for (c) and (d). Filled circles and open triangles correspond 
to the ground state and excited states, respectively. 
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FIG. 6: Spin gap (solid line) and binding energy (broken line) 
as a function of n (= 1 — Nh/L): 7 — 0.6 and Ji/ti — 0.25. 
The dotted lines with filled circles are the results of As for 
L = 24, 32, 40, 48, 56, 64, 72, 80 from top to bottom. 



IV. SPIN EXCITATIONS IN A METALLIC 
PHASE 

We now discuss spin excitations in a hole-doped system 
near half filling. We start with low-energy excitations in 
a system with two holes. As mentioned above, the sys- 
tem is driven to a metallic phase immediately upon hole 
doping for the dimer phase (0 < 7 < 0.71). On the other 
hand, for the plaquette phase (7 > 0.71), doped holes 
may be in a localized state, as a result of a first-order 



transition. Namely, holes are accommodated in a degen- 
erate dispersionless level up to a critical hole density. In 
the following, we thus pay our attention to the former 
region, < 7 < 0.7, which is more interesting as far as 
the system with a small number of holes is concerned. 

In order to discuss excitations in a doped case, we con- 
sider two quantities, i.e. the spin gap As{Nfi) and the 
binding energy As(-/Vh) respectively defined for (even 
number) holes on 2 x L lattice sites as, 

As{Nh) = Eo{Nh, S' = l)- Eo{Nh, = 0), (3) 
AeiNh) = 2Eo{Nh-l,S^ ^l/2)-Eo{Nh,S' ^0) 
-Eo{Nh-2,S^ = 0). (4) 

In Fig. ^, we show the spin gap A5 and the binding 
energy Ab computed for two doped holes {Nh ~ 2). The 
inset shows a finite-size scaling of A5 and As for 7 = 0.6. 
It should be noticed that the spin gap and the binding 
energy monotonically increase as 7 increases from 7 ~ 
0.35. It is seen that the spin gap A5 just coincides with 
the binding energy Ab, implying that a pair-breaking 
excitation in a metallic phase may be described by two 
independent quasiparticles. It is to be noted that the 
above spin gap in a doped system is totally different from 
the spin gap defined at half filling, A5 = Ji, as we will 
see below. 

Shown in Fig. || is a spatial distribution of the hole 



density (nf ) {n'- 



1 



and the spin density (Sf) for 



two holes on a L = 64 system. Here filled circles and open 
triangles show the results for the singlet ground state and 
for a triplet excited state, respectively. For 7 = 0.6, as 
shown by filled circles in Fig. H(c), a hole-pair exists in 
the ground state, implying the existence of an attractive 
interaction between two holes in the ground state. For 
a triplet excitation, the bound state is broken and two 
quasiparticles appear, as seen from open triangles in Fig. 
^(c) and (d). On the other hand, the situation is some- 
what different for 7 — 0.3, where the repulsive interac- 
tion seems to exist even in the ground state, as shown by 
filled circles in Fig. ||(a) and (b). We have checked that 
a repulsive interaction in the latter case is mainly caused 
by finite-size effects. Therefore, it is expected that such 
a repulsive interaction becomes small for large L, and 
in the thermodynamic limit, doped holes may behave as 
two independent quasiparticles in the case of 7 = 0.3. 

We now turn to the case with a finite density of holes. 
In Fig. H, we show the spin gap as a function of the elec- 
tron density n for 7 = 0.6. The dotted lines with filled 
circles are the results of A5 for finite-size systems. It is 
seen that the value of the spin gap dramatically changes 
between half-filling and two-hole case, and then smoothly 
changes in the region Nh > 2. We have performed a 
finite-size scaling analysis of the spin gap from the data 
shown by dotted lines, obtaining the curve shown by the 
solid line. Note that the spin gap has a discontinuity 
at half filling, implying that the spin gap in a metal- 
lic phase has a different origin from the half-filled case. 
Also shown by the broken line is the binding energy A^ 
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computed by a similar finite-size scaling analysis. Since 
A5 is nearly equal to A^, we can say that excited states 
may be described by two independent quasiparticles, as 
we have already seen in a system with two holes. 

In order to see the nature of the attractive interaction 
between holes more clearly, we compute short-range hole- 
hole correlation functions Ha in a bulk limit, which are 
defined in the rung (_L) and chain (||) directions. 



h\2 



In Fig. ^ we show Ha as a function of n for 7 — 0.3 and 
0.6. For 7 — 0.6, there is the region with both of H± 
and being positive, where an effective interaction be- 
tween holes is attractive in both directions. On the other 
hand, for 7 = 0.3, both of H± and H\\ are negative (i.e. 
repulsive interaction) in the whole range shown in the 
figure. By comparing these results with Fig. ^, we see 
that the region, where holes attract each other, roughly 
corresponds to the spin-gap metallic region. This means 
that the mechanism of spin-gap formation in a metal- 
lic phase, which is different from that in the insulating 
phase, is closely related to the attractive force between 
holes. 



H -0.01 




FIG. 7: Short-range hole-hole correlations Ha for 7 = 0.3, 0.6 
and Ji/ti = 0.25. 



We have found that the transition from the insulat- 
ing plaquette phase to a metallic phase is of first-order, 
being contrasted to the dimer phase showing a second- 
order transition. Furthermore, a spin-gap metallic phase, 
which is indicated by the shaded region in Fig. ^, per- 
sists in the vicinity of half-filling, which is indeed caused 
by strong frustration. For the value Ji/ii — 0.25 shown 
in Fig. 1^, we have found that a spin-gap metallic phase 
appears in the region of 7 > 0.35(= 75). 



When the electron density is further decreased, we en- 
counter a CDW insulating phase at quarter filling, char- 
acterized by the plateau in n-/i curve. The region of a 
CDW phase becomes narrower as 7 increases, and finally 
vanishes beyond a certain critical value. As is seen from 
Fig. ||(b) and Fig. ||(b), a metallic phase continuously 
spreads in the whole region of 7, where short-range spin 
correlations along the rung, S_l, always take negative val- 
ues. This implies that dimer correlations are dominant 
for the ground state in a metallic phase, as is usually 
the case for a doped antiferromagnetic ladder. However, 
it should be noticed that in a metallic phase close to 
the plaquette-singlet insulator, a competition between 
the dimer- and plaquette- singlet states may give rise 
to dual properties in physical quantities. 



We have thus far shown the calculated results by choos- 
ing a specific ratio of hopping and exchange. It is neces- 
sary to check what is changed when we take other choices 
of the parameters. For instance, we show the spin gap 
of a system with two holes in Fig. ^ when the ratio of 
Ji/ii is changed. As Ji/ti decreases, the critical point, 
7s, where the spin gap starts to develop, is shifted to a 
larger value. Note, however, that 7s 0.5 in the limit of 
Ji ^ ti. Namely, even if Ji/ti is small, a spin-gap metal- 
lic phase appears at least in the region of0.5<7<0.71. 
This feature is diffepini,-fcaaiJJ^ case of the t-J ladder 
without frustrationjHjOEj'E3ot3 which has a spin gap 
only for relatively large Ji/ti. 



V. PHASE DIAGRAM 



Based on the above results, we now obtain the ground- 
state phase (diagram, which is shown in Fig. ^j. As found 
by Gelfandjtj the ground state at half filling is either in 
the dimer-singlet (or plaquette-singlet) phase with spin 
gap at half filling for 7 < 0.71 (7 > 0.71). These two 
phases are separated via a first-order phase transition 
when 7 is changed. 
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FIG. 8: Ground-state phase diagram for J\/ti — 0.25: (a) 
dimer-singlet phase (n = 1), (b) plaquette-singlet phase (n — 
1), (c) CDW phase (n = 1/2), (d) empty phase (n = 0), (e) 
metaUic phase, and shaded area indicates a spin-gap metallic 
phase. 




FIG. 9: Spin gap with two holes as a function of 7 for 
Ji/ti = 0.25,0.10,0.05. Other parameters are chosen as in 
Fig. §. 



to a metallic phase. However, when the system is in 
the plaquette phase, a first-order transition occurs upon 
hole doping. The latter first-order transition occurs as a 
consequence of a strong competition between the dimer- 
and plaquette-singlet states. This is different from ordi- 
nary metal- insulator transitions for a frustrated J chain 
with_next-nearest neighbor exchange couplings studied so 
far.N 

The first-order phase transition found here is closely re- 
lated to unique geometry of the linked-tetrahedra lattice, 
which causes a dispersionless (spatially isolated) mode, as 
mentioned in the beginning of the paper. It is to be noted 
that such a dispersionless mode caused by special geom- 
etry exists for the 3D pyrochlore lattice and also for the 
2D orthogonal-dimer lattice. In this sense, the present re- 
sults may capture some of characteristic properties com- 
mon in this class of fully frustrated systems with doped 
holes. For instance, a pyrochlore antiferromagnetic sys- 
tem, indeed shows a competition between the dimer and 
plaquette phases, being similar to the present results. eI 
Therefore, we think that unusual properties may ap- 
pear in a metallic phase of the doped pyrochlore system 
around the insulating phase. This is consistjsnt with the 
results deduced by preservative analyses. Dtil We also ex- 
pect that such peculiar properties may emerge even in 2D 
orthogonal-dimer systems when holes are doped into an 
insulating phase around the boundary between the dimer 
and plaquette phases.Q Detailed theoretical studies on 
hole-doped pyrochlore systems as well as 2D orthogonal- 
dimer systems are left as an interesting open problem. 



VI. SUMMARY AND DISCUSSIONS 
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